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As very high energy (> 100 GeV) gamma rays travel over cosmological distances, their flux is 
attenuated through interactions with the extragalactic background light. Observations of distant 
gamma ray sources at energies between ~2Q0 GeV and a few TeV by ground-based gamma-ray tele- 
scopes such as HESS, however, have motivated the possibility that the universe is more transparent 
to very high energy photons than had been anticipated. One proposed explanation for this is the 
existence of a^xion-like-particles (ALPs) which gamma rays can efhciently oscillate into, enabling 
them to travel cosmological distances without attenuation. In this article, we use a state-of-the-art 
model for the extragalactic background light (which is somewhat lower at ~ /im wavelengths than 
in previous models) and data from the Fermi Gamma Ray Space Telescope to calculate the spec- 
tra at 1-100 GeV of two gamma-ray sources, lESllOl-232 at redshift z = 0.186 and H2356-309 at 
z = 0.165, in conjunction with the measurements of ground-based telescopes, to test the ALP hy- 
pothesis. We find that these observations can be well-fit by an intrinsic power-law source spectrum 
with indices of -1.72 and -2.1 for lESllOl-232 and H2356-309, respectively, and that no ALPs or 
other exotic physics is necessary to explain the observed degree of attenuation. While this does not 
exclude the possibility that ALPs are involved in the cosmological propagation of gamma rays, it 
does reduce the motivation for such new physics. 

PACS numbers: 95.85.Pw, 98.70.Vc 98.70.Rz 14.80.Mz 
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I. INTRODUCTION 

Very high energy gamma rays {E^ > 100 GcV) scat- 
tering with photons in the infrared to ultraviolet range 
can exceed the threshold for electron-positron pair pro- 
duction. As a consequence, gamma rays are predicted to 
experience significant attenuation over cosmological dis- 
tances, leading to the suppression of the gamma-ray flux 
from high redshift sources such as active galactic nuclei 
and gamma-ray bursts [1]. 

Beginning several years ago, however, ground-based at- 
mospheric Cerenkov telescopes (HESS, MAGIC, VERI- 
TAS) began to report measurements of distant gamma- 
ray sources which did not appear to exhibit as much at- 
tenuation as had been generally anticipated [2, 3]. In 
particular, the HESS collaboration reported spectra from 
two relatively high redshift sources, the active galactic 
nuclei (AGN) H2356-309 at z = 0.186 and lESllOl-232 
at z = 0.165, that appeared to exhibit relatively little at- 
tenuation. These observations indicate that either these 
particular sources are injecting a very hard spectrum of 
gamma rays {dN^/dE^ oc E~^ with F ~ —0.6 to 0.2, in 
contrast to F '^ 2 as predicted from second order Fermi 
acceleration), or that the universe is considerably more 
transparent to very high energy photons than had been 
expected [4, 5]. If the latter is the case, it could indicate 
a low density of extragalactic background light (EBL) 
or, alternatively, the existence of exotic physics which 



enables very high energy photons to traverse over cosmo- 
logical distances without experiencing attenuation. 

One possibility that has received some attention is that 
very high energy photons may be oscillating into axion- 
like particles (ALPs), which do not experience significant 
attenuation, only to be reconverted into photons before 
reaching Earth. In Ref. [6] (and later Ref. [7]) it was 
proposed that the magnetic fields in and around the very 
high energy gamma-ray sources could facilitate the ef- 
ficient conversion of very high energy gamma rays into 
ALPs. This effect would be expected to become signifi- 
cant at the GeV- TeV energy range currently being stud- 
ied by satellite and ground-based gamma-ray telescopes. 
Furthermore, it is plausible that a significant fraction of 
such particles could be reconverted into photons through 
interactions with the magnetic field of the Milky Way [8] . 
Alternatively, it was suggested in Ref. [9] that photon- 
ALP conversions could take place through interactions 
with the (nano-gauss scale) intergalactic magnetic field. 
Another explanation consistent with both the EBL cal- 
culations and AGN models involves ultra high energy 
gamma rays of energies below 50 EeV crossing cosmo- 
logical distances and interacting with the EBL close to 
Earth, thus generating secondary photons that are ob- 
served by gamma-ray telescopes [10, 11]. 

The measurements of the spectra of cosmologically dis- 
tant gamma-rays sources by ground-based atmospheric 
Cerenkov telescopes have thus far been unable to distin- 



guish between the various possible explanations for the 
lack of observed attenuation in these sources. Studies 
of these sources' spectra over a wider range of energies, 
however, could be used to constrain their initial (unat- 
tenuated) spectral shapes, and to determine with greater 
precision the degree to which very high energy gamma 
rays are attenuated by the EBL. In particular, because 
the EBL has a negligible density at wavelengths below 
0.1 /zm, there is no significant attenuation of gamma rays 
with energies below the corresponding threshold energy 
of ^20 GeV, and only modest attenuation below ^100 
GeV. 

The most significant source of uncertainty in deter- 
mining the attenuation of gamma rays is the EBL spec- 
trum. The two complementary methods of estimating 
the EBL spectrum are the integrated galaxy counts, de- 
pendent on star formation history, and the absolute mea- 
surements, prone to uncertainties due to foreground con- 
tamination. As it stands the model proposed by Frances- 
chini et al. [12] provides the minimal spectrum for the 
first peak in the optical to near infrared (NIR) and cor- 
responds to minimal values of the calculated total EBL 
[13]. In this article, we use publicly available data from 
the Large Area Telescope (LAT) on the Fermi Gamma 
Ray Space Telescope to study the spectrum of two active 
galactic nuclei (lESllOl-232 and 112356-309) between 1 
GeV and 100 GeV. When combined with the higher en- 
ergy spectral measurements of ground-based gamma-ray 
telescopes, a more complete picture is obtained, cover- 
ing more than three decades of energy, including a range 
for which there is expected to be little or no attenuation 
due to electron-positron pair production. We find that 
these measurements imply a relatively low density of the 
EBL, but that the observed spectra can be accounted for 
without the presence of ALPs or other exotic physics. 



II. THE ATTENUATION OF COSMOLOGICAL 
GAMMA-RAY SOURCES 

As very high energy gamma rays propagate through 
the universe, their interactions with the EBL and the 
resulting production of electron-positron pairs cause their 
spectrum to become attenuated by a factor of e^'^(-^""^°\ 
where the optical depth, t{E~^, zq), is given by 






with 
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do; 



(1) 



a{E^,u;,z) = -(l-/x)cr^^(£'^,a;,^) dpi . 

(2) 
Here, zq is the redshift of the source, H{z) is the rate of 
Hubble expansion, E^ and w are the (appropriately red- 
shifted) source and background photon energies, respec- 
tively, fi is the cosine of the angle between the incoming 



and target photon, a^^y is the cross section for electron- 
positron pair production, and driEBL/dw denotes the 
spectrum and density of the EBL [14]. 

For a gamma ray of energy E^, the pair production 
cross section is largest for background photons of energy 



c^„,ax ^ 2(TOec') V^7 - 0-5 (^^) eV, (3) 

or equivalently 



1.24 ^m X 



E'^ 



ITeV 



(4) 



From this expression, we see that gamma rays with 
Ej — 100 GeV (approximately the highest energies mea- 
sured by Fermi) interact most efficiently with the far-UV 
component of the EBL. Lower energy gamma rays, con- 
stituting the bulk of the Fermi events, experience rela- 
tively little attenuation through pair production, since 
the EBL is negligible at far-UV and shorter wavelengths. 
In contrast, the HESS measurements of gamma rays at ~ 
2 TeV directly probe the EBL in the near-IR and longer 
wavelengths. 

The EBL consists of the photons emitted by extra- 
galactic sources, both resolved and unresolved, over all 
cosmic epochs. Its density and spectrum are difficult 
to measure directly due to the presence of bright at- 
mospheric and galactic foregrounds, including Zodia- 
cal scattered light, interplanetary dust emission. Galac- 
tic starlight, and high Galactic latitude cirrus emission 
(Galactic stellar light reflected by high latitude dust). In 
Fig. 1, we show three examples of EBL models, as de- 
scribed in Refs. [2, 12, 15]. 

The distribution of extragalactic optical and infrared 
background photons used in this paper is that from 
Franceschini et al. [12]. Their analysis of the EBL 
uses survey data in the optical, near-IR, and sub- 
millimeter from ground-based observatories as well as 
multi-wavelength information from Hubble Space Tele- 
scope (HST), Infrared Space Observatory (ISO), and 
Spitzer Space Telescope, while also utilizing direct mea- 
surements of and upper limits on the EBL by COBE 
to place additional constraints on their estimate. The 
Franceschini et al. model of the optical background is 
based on HST deep galaxy count estimates. To avoid 
the large near-IR excess found in some analyses based 
on the COBE measurements at 1.2 and 3.5 ^m, and 
the accompanying strong discontinuity in the EBL at 
~ 1/xm, they use a model of the near-IR background de- 
rived from direct integration of the ultra-deep HST and 
Spitzer counts. 

Referring to Fig. 1, we point out that the model of 
Franceschini et al. [12] has a considerably lower estimate 
for the EBL's density in the near-IR, and a somewhat 
higher estimate in the far-IR, relative to the other mod- 
els shown. Franceschini et al. [12] constrain the near-IR 
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FIG. 1. Representative models of the extragalactic back- 
ground light (EBL) compared to various measurements. The 
solid black line is the model of Franceschini et al. [12] , which 
we adopt in our calculations. The dashed red line denotes the 
model used by the HESS collaboration in Ref. [2] . The dotted 
blue line is the model of Primack et al., using the Kennicutt 
et al. stellar initial mass function [15, 16]. The data are as 
described in Ref. [12]. The three rightmost data points in 
the far-IR are from the analysis of DIRBE data by Madau et 
al. (black squares) [17] and the re-analysis of that data by 
Lagache et al. (red squares) [18]. The two lower limits given 
at 70 and 160 ^m (cyan filled circles) come from a stacking 
analysis of Spitzer data described in Ref. [19]. The mid-IR 
data points at 15 and 24 /xm are the resolved fraction of the 
EBL by the deep ISO surveys IGTES [20, 21]. The near-IR 
data point at 8 ^m (open triangle) is from Ref. [12], while 
the three lower wavelength points (open triangles) are from 
Ref. [22], using the integration of the IRAC galaxy counts. 
The far-UV/optical/near-IR points from 0.25 ^m to 2.5 /im 
(open squares) are from an analysis of Ref. [17] based on deep 
galaxy counts. 



at 8 /im using Spitzer IRAC observations of the faint 
number counts [23, 24] and the local galaxy luminosity 
function [25] . These data significantly constrain the EBL 
at the junction between the stellar photospheric emission 
in the optical/near-IR and the far-IR dust emission peak. 
The most relevant yet controversial region of EBL for 
the attenuation of gamma rays of energies of 100-500 GeV 
is the range from of optical through near IR. The anal- 
ysis by Bernstein et al. [26] obtained background values 
a factor of ^ 3 greater than those estimated from direct 
HST counts. Their approach was based on subtraction 
the contributions of the Zodiacal light and the diffuse 
galactic light from the total signal measured by HST. 
Such a discrepancy would imply that either there exists 
a truly diffuse background of unknown origin or that the 
HST missed a large part of the light from the faint galax- 



FIG. 2. The attenuation factor for gamma rays due to 
electron-positron pair production for sources at redshifts of 
z = 0.1, z = 0.2, z = 0.3, z = 0.4, and z = 0.5, from right- 
to- left. The dashed and solid curves denote the attenuation 
for the EBL models of Franceschini et al. [12] and Primack et 
al. [15, 16], respectively (see Fig. 1). 



ies. It was later demonstrated that the original analysis 
strongly underestimated the systematic errors and that 
the background estimates have to be considered as upper 
limits [27]. 

Throughout our analysis, we assume that the shape 
and normalization of the background radiation docs not 
change except for rcdshifting between the source and the 
observer. As suggested by the models of Ref. [28], this 
is a reasonable assumption for z < 0.4. In Fig. 2, we 
plot the attenuation resulting from this EBL model, as 
a function of observed energy, and for sources at various 
redshifts. The energy dependence of the optical depth 
leads to modifications in the source spectrum as mea- 
sured locally. In general, even a weak dependence on en- 
ergy of the optical depth can lead to large changes in the 
source spectrum because of the exponential dependence 
of the attenuation. Since the optical depth is generally 
larger for photons with higher energies, the effect of at- 
tenuation is to make the locally observed spectrum softer 
than the source spectrum, as well as to reduce the overall 
intensity. 



III. THE MODIFICATION OF GAMMA-RAY 
SPECTRA BY AXION-LIKE PARTICLES 

As discussed in Refs. [6, 7], the probability of a pho- 
ton converting into an ALP in a magnetic field is related 
to the size of the domain of the magnetic field and the 



strength of the magnetic field component along the polar- 
ization vector of the photon. Given the typical sizes of as- 
trophysical accelerators and the magnetic field strengths 
present in them, significant photon-to-ALP conversion 
can occur in or near very high energy gamma-ray sources 
over a large range of allowed ALP parameter space. For 
an unpolarized photon source, one expects 1/3 of the 
original photons to be converted into ALPs at the source 
above a critical energy, £. That is, for energies E^ ^ £ 
there could potentially be an ALP flux from very high 
energy gamma-ray sources as large as ^^50% of the resid- 
ual gamma-ray flux exiting the source. The critical en- 
ergy, £, can naturally fall in the gamma-ray band. The 
photon- ALP mixture then propagates toward the Milky 
Way without significant further oscillations, since in the 
absence of an external electromagnetic field, the ALPs 
and the photons effectively decouple. During propaga- 
tion over cosmological distances, the spectrum of pho- 
tons is depleted via pair production. As a result, upon 
reaching the Milky Way, the very high energy flux will 
be dominated by ALPs. 

The reconversion probability of ALPs into photons in 
the the Galactic Magnetic Field depends on a number 
of factors, including the ALP-photon coupling, the ALP 
mass, the photon energy, and the strength and geometry 
of the magnetic field. The CERN Axion Solar Telescope 
(CAST) experiment [29, 30] provides a direct bound on 
the ALP-photon coupling of ga-y < 8.8 x lO^^GeV^^ 
for rua < 0.02 eV. For ultra-light ALPs {rria < 10"" eV), 
the absence of gamma rays from SN 1987A yields a limit 
of Sfa7 < 1 X 10-"GeV-i [31] or even gaj < 0.3 x 
10-"GeV"i [32]. However, in the range of 10"" eV < 
TUa ^ 10^^ eV, the CAST bound is the most general and 
stringent. The optimal range of parameters for the re- 
conversion mechanism is 10"^" < nia < 10~^, over which 
couplings as large as ^q^ w 8 x 10~^^GeV~^ are consis- 
tent with present bounds. Although the characteristics 
of the Galactic Magnetic Field are not very well known, 
especially in the directions toward the Galactic Center 
or away from the Galactic Plane, reasonable models of 
the Galactic Magnetic Field predict regions of the sky in 
which the reconversion probability can be 20% or larger 
for an ALP-photon coupling of 5a = 5 x 10~^^GeV~^ [8]. 
Moreover, Rcf. [8] found that, in these models, a sig- 
nificant fraction of the sky corresponds to a probability 
larger than 10%. Thus, appreciable reconversion proba- 
bilities are possible, though difficult to predict given the 
scarcity of our knowledge regarding the structure of the 
Galactic Magnetic Field. 

We note that for the most viable range of parame- 
ters for which the ALP mechanism is efficient (ga ~ 
5 x lO^^^GeV"-^, rUa ^ 10"^ eV and microgauss-scale 
fields, for example), the critical energy, £, above which 
the probability of a photon converting into an ALP is of 
the order unity, falls in the sub-GeV range. Therefore, 
the observed gamma-ray spectrum will not contain pe- 
culiar features resulting from the onset of photon- ALP 
oscillations. Thus, over the whole range of energies ob- 



served by Fermi and atmospheric Cherenkov telescopes, 
the ALP flux generated at the source via oscillation will 
have the same spectrum as the photons at the source. 

In Fig. 3, we show the effects of ALPs on the gamma- 
ray spectrum at Earth for the case of a source spectrum 
of dN/dE oc E~^ with F = 2 after propagating from 
z = 0.2 (top) and z = 0.5 (bottom), assuming oscillation 
into ALPs at the source with a probability of 1/3 and 
reconversion in the Galactic Magnetic Field with a prob- 
ability of 0.1. For comparison, we also include the results 
for a standard scenario with no ALPs. Clearly, in the 
presence of ALP-photon mixing, the quasi-exponential 
cutoff is interrupted, leading the spectrum to plateau at 
high energies. 



IV. THE MEASURED SOURCE SPECTRA 

The optical depth of the universe to very high en- 
ergy photons increases with redshift, so the spectra of 
very distant sources experience greater attenuation. Any 
mechanism mitigating the attenuation is therefore most 
apparent in the spectra of high redshift sources. In 
our analysis, we consider two relatively high redshift 
gamma-ray sources which have each been observed by 
both HESS [2] and Fermi [33], lESllOl-232 and H2356- 
309. These sources arc both blazars, and were flrst 
detected at TeV energies by HESS in 2006. lESllOl- 
232 is located at a distance oi z = 0.186 at {l,b) = 
(273.19°, 33.08°). H2356-309 is located at a distance of 
z = 0.165 at {l,b) = (12.86°, -78.04°). HESS has re- 
ported data on the spectra of these sources in the energy 
range - 200 GeV < £■ < 3 TeV [2]. Fermi has not de- 
tected variability in either of these sources [33] . 

The Fermi-LAT team made their first year gamma-ray 
event data publicly available at the end of August 2009, 
and this data is now updated daily on their website [34] . 
Additionally, Fermi has supplied the Fermi Science Tools 
for the analysis of LAT data [35]. In our analysis we 
use the LAT event data collected over the period from 
August 4, 2008 to January 11, 2010. 

The LAT point spread function (PSF) for front- 
converting events^ is a decreasing function of photon 
energy with a 95% containment radius of ~1.2° at 1 
GeV and less than 0.25° above 10 GeV. Back-converting 
events have a PSF of --2° at 1 GeV and less than 0.5° 
above 10 GeV. We note that an analysis by Ref. [36] 
suggests that the width of the PSF for back-converting 
photons with energies above ~10 GeV may be underes- 
timated in the Monte Carlo simulations used to gener- 
ate the LAT response functions. To prevent contami- 
nation of the spectra by other nearby sources, we limit 



^ The front, or "thin", section of the LAT tracker has twelve 
tungsten converters, each of which is 0.035% radiation lengths 
thick, while the back, or "thick", section has four 0.18% radia- 
tion length tungsten converters. 
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FIG. 3. The gamma-ray spectrum at Earth from a source 
with an injected spectrum of dN/dE oc E~^, after propaga- 
tion over a distance of 2 = 0.2 (top) and z = 0.5 (bottom). 
Results are shown with and without the effects of photon- ALP 
oscillations. The ALP-photon mixing mitigates the impact of 
absorption via pair production and leads to a plateau in the 
spectrum at high energies. We have calculated the effects of 
ALP-photon mixing assuming a source conversion probability 
of 1/3 and a Milky Way reconversion probability of 0.1. The 
vertical axis is in arbitrary units. 



the angular size of our analysis region according to a 
radius equal to half of the angular distance to the near- 
est known source. For lESllOl-232 and H2356-309, the 
closest known sources arc 239.151' and 125.244' away, 
respectively. The energy dependence of the PSF thus 
determines the minimum energy in the analysis of the 
spectrum. For consistency we calculate the spectra for 
both AGN over the same energy range, 1-100 GeV. 
To determine the spectra of the AGN, we calculate 
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FIG. 4. The spectra of gamma-ray sources lESllOl-232 at 
z = 0.186 (upper) and H2356-309 at z = 0.165 (lower) as 
determined from Fermi data {black circles) and as measured 
by HESS {black squares). The curves shown are the calcu- 
lated best fits to the data assuming absorption of gamma 
rays via electron-positron pair production interactions with 
the EBL. The intrinsic source spectra are taken to be of the 
form dN/dE oc E~^ . The best-fit values for the power law 
indices are F = 1.72 for lESllOl-232 and F = 2.1 for H2356- 
309. The EBL model used is that of Franceschini et al. [12]. 
See text for more details. 



the total flux and subtract off the background. We esti- 
mate the background spectrum using two methods, which 
we refer to as the "Fermi" method and the "annulus" 
method. For the Fermi method, we assume that the back- 
ground flux is the sum of the diffuse galactic emission [37] 
and the isotropic emission [38] (which is also referred to 



as the extragalactic diffuse emission) as measured by the 
LAT team. The diffuse flux has its source in cosmic rays 
interacting with the interstellar gas and interstellar radi- 
ation field, while the isotropic flux is composed mainly of 
gamma rays from unresolved extragalactic sources. For 
the annulus method, we calculate the total flux in an 
annulus of width 0.25°, centered around the source, and 
with an inner radius that is 10% larger than the PSF at 
a given energy. 

These two methods for calculating the background flux 
in the region of lESllOl-232 give results that agree very 
well. However, for H2356-309, the Fermi and the annu- 
lus backgrounds do not agree; the annulus background 
is more than an order of magnitude greater than the 
Fermi background at ~2 GeV, and exceeds the Fermi 
background by even larger factors at higher energies. 
We suspect that this can be explained by the exis- 
tence of a currently unresolved point source in the re- 
gion and/or by contamination from the nearby source 
IFGL J2350. 1-3005 due to the underestimation of the 
PSF suggested by Ref. [36], as mentioned earlier in this 
section. The Fermi spectra shown in the upper frame 
of Fig. 4 (lESllOl-232) was calculated using the Fermi 
background method, while the results in the lower frame 
(H2356-309) have made use of the annulus method. The 
error bars shown for the spectra in Fig. 4 were calculated 
assuming statistical uncertainties in the total flux and 
systematic-|-statistical uncertainties (though systematics 
dominate) in the isotropic flux. 



law index of F w 0.6 in the absence of ALPs. The large 
discrepancy between the values of F calculated by Ref. [8] 
(-0.2 and 0.6) and the values found in this study (1.72 and 
2.1) is largely due to the difference in the EBL models 
used. For gamma rays with energies in the HESS range, 
the dominant contributors to pair production are EBL 
photons with wavelengths of 3.75 /xm > A > 0.25 /im. 
The density of the EBL in the model used in our cal- 
culations [12] is a factor of ^2 — 3 times smaller than 
that used in Ref. [8] (see Fig. 1). Thus, the predicted at- 
tenuation is much less and the required intrinsic source 
spectrum much softer than that found in Ref. [8] . More- 
over, the Fermi data at energies of '--'l-lOO GeV confirms 
this conclusion by constraining the spectrum in a energy 
range largely unaffected by attenuation. 

A similar analysis was done in Ref. [41], where 
lESllOl-232 and H2356-309 were analyzed, among other 
gamma-ray sources, using HESS but not Fermi data. 
They used two EBL models of Ref. [42] corresponding 
to different galaxy evolution models, a baseline model 
(B) and a fast evolution model (FE). The calculated in- 
trinsic power-law index of lESllOl-232 in their analysis 
was found to be F w 1.0 using model B and F « 1.5 
using model FE. For H2356-309 the intrinsic power-law 
index was found to be F « 1.5 and F p» 1.9, for models B 
and FE, correspondingly. These values are much closer 
than those of Ref. [8] to our present results, which can 
explained by the similarities of the EBL models. 



RESULTS 



In Fig. 4, we compare the spectra of lESllOl-232 and 
H2356-309 as observed by Fermi and HESS to the shape 
predicted after propagation, assuming an intrinsic source 
spectrum of the form dN^/dE oc E~^ and using the EBL 
model of Ref. [12] (see Fig. 1). We find that in the case of 
lESllOl-232, the combination of Fermi and HESS data 
is best-fit for a value of F = 1.72 (x^/dof « 0.43). For 
H2356-309, a value of F = 2.1 yields an acceptable fit 
to the spectrum above ^3 GeV {x^/dof « 0.61). One 
should keep in mind that the injected spectrum from 
astrophysical sources need not follow a simple power- 
law [39, 40], although power-law forms do yield good fits 
to the observed spectra in the cases considered here. On 
this basis, we conclude that the observed spectra of these 
sources do not contain compelling indications of axion- 
like particles. 

A few comments are in order at this time. In Ref. [8], a 
similar calculation, for the same two gamma-ray sources 
analyzed here, lESllOl-232 and H2356-309, was per- 
formed using HESS (but not Fermi) data and the EBL 
model of Ref. [15]. In their analysis, they found that, un- 
der standard astrophysical assumptions (ie. no ALPS), 
the intrinsic source spectrum of lESl 101-232 was re- 
quired to have a very hard power-law index of F ~ —0.2. 
Similarly, they found that H2356-309 required a power 



VI. SUMMARY AND CONCLUSIONS 

In this article, we have studied the spectra of two 
gamma-ray sources, lESllOl-232 at redshift z ~ 0.186 
and H2356-309 a,t z — 0.165, based on the observations 
of the the Fermi Gamma Ray Space Telescope and the 
ground-based telescope HESS. While the observations of 
HESS at energies above a few hundred GeV indicate a 
somewhat surprising degree of the transparency of the 
universe to very high energy photons, the reason for this 
had not until now been clear. By including in our analysis 
recent data at lower energies (<100 GeV) from Fermi, we 
find that the spectra of lESllOl-232 and H2356-309 are 
consistent with an intrinsic source spectrum with power 
law indices of approximately -1.7 and -2.1, respectively. 
The lack of attenuation observed by HESS can be ac- 
counted for if the density of the extragalactic background 
light (EBL) is somewhat low in the UV-optical-near-IR 
range (such as that in the model of Franceschini et al). 

Having emphasized the importance of the EBL model 
we conclude that no exotic physics, such as oscillations 
between photons and axion-likc particles (ALPs), is re- 
quired to explain the relative lack of attenuation observed 
by HESS. 
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